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ABSTRACT 



Aims. This work reports radiative transition rates and electron impact excitation rate coefficients for levels of the n= 
3, 4, 5, 6, 7, 8 configurations of Call. 

Methods. The radiative data were computed using the Thomas-Fermi-Dirac central potential method in the frozen 
core approximation and includes the polarization interaction between the valence electron and the core using a model 
potential. This method allows for configuration interactions (CI) and relativistic effects in the Breit-Pauli formalism. 
Collision strengths in LS-coupling were calculated in the close coupling approximation with the R-matrix method. Then, 
fine structure collision strengths were obtained by means of the intermediate-coupling frame transformation (ICFT) 
method which accounts for spin-orbit coupling effects. 

Results. We present extensive comparisons with the most recent calculations and measurements for Call as well as a 
comparison between the core polarization results and the "unpolarized" values. We find that core polarization affects the 
computed lifetimes by up to 20%. Our results are in very close agreement with recent measurements for the lifetimes of 
metastable levels. The present collision strengths were integrated over a Maxwellian distribution of electron energies and 
the resulting effective collision strengths are given for a wide range of temperatures. Our effective collision strengths for 
the resonance transitions are within ~11% from previous values derived from experimental measurements, but disagree 
with latter computations using the distorted wave approximation. 

Key words, atomic data - atomic processes - line: formation - stars: Eta Carinae - Active Galactic Nuclei 



' 1 1. Introduction 



Call plays a prominent role in astrophysics. The so-called 
H and K lines of this ion ar e important probes of so - 
lar and stellar chromospheres (jRauscher fe Marcvl [2006) . 
In the red spectra of Active Galactic Nuclei (AGN) the 
infrared triplet of Can in emission has b een used t o in- 
vestigate the correlations with optical Fe II (|Jolvlll989f ) and 
thei r implications on the phys ical conditions of the emitting 
gas (jFerland fc Perssonlll989l ) . [Ca n] optical emission lines 
together with the infra red [Fell] are often used as probe of 
dust content of AGNs (jShields et al.lll999t ). 

Can has been addressed by numerous theoretical and 
experimental groups. The lifetimes r for the 4p 2 P° and 
3d 2 D levels have been measured wi th high precision 
(|Jin fc Churchlll993HKreuter et"al"1 l2005). Various theoret- 
ical methods have been used in trying to match the se ex- 
perim ental values. The most recent calculations of ILiawl 
(1995) using the Brueckner approximation with third-order 
correction agree within ~1% with the experimental life- 



* The atomic data 
ergy levels, A- values, 
is available 
mous ftp 



from this work, including en- 
and effective collision strengths, 
in electronic form at the CDS via anony- 
to cdsarc.u-strasbg.fr (130.79.128.5) or via 



http:/ /cdsweb. u-strasbg.fr/cgi-bin/qcat?J/A-|-A/ 



time for the 4p 2 P° levels but are ~11% to o small for the 
3d 2 D metastable levels. The calculations of iKreuter et al.l 
(pool , using a relativistic all-order method which sums in- 
finite sets of many-body perturbation theory terms, agree 
within ~0.3% with the experimental lifetimes of the 3d 2 D 
metastable levels of Ca ll, but offer no data for other levels. 
iGuet fc Johnson! (|l99lh computed lifetimes using relativis- 
tic many-body perturbation theory that agree within ^2% 
with experimental values for the Ap 2 P° levels and within 
^6% for the metastable states. In the present work we use 
the Thomas-Fermi-Dirac central potential with core polar- 
ization interaction to provide a complete set of accurate 
A-values for allowed and forbidden transitions to be used 
in modeling Can spectra. 

Various calculations of collision strengths have been 
performed for the re s onance transition s in Call (see 
I Zatsarinnv et all Il991t IChidichimol 119811; iKennedv et alj 
Il978t ISaraphl I1970D . lOsterbrock fc Wallace! (|1977l) de- 
rived effective collision strengths from experimental cross 
sections of the resonanc e K and H lin e s of Can at 

3934 and 3968 A by iTavlor fc Dunr] (|1973D . Later, 

IZapesochnvi et all (|1975f) published cross sections for ex- 
citing 5s and 4d levels from the ground state, which are 
importa nt to estimate the c ontribution to the 4p level by 
cascade. Mi trov et al.1 (|1988l ) presented a detailed study of 
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the electron- impact excitation of the (4s-4p) transitions us- 
ing the close-coupling approx imation includ i ng a p olariza- 
tion potential. More recently [Burgess et al.1 (|l995t ) used a 
non-exchange distorted wave approximation including the 
lowest 7 Can terms. 

The IRON Project is an international enterprise 
devoted to the computati on of accurate atomi c data for 
the iron group elements (jHummer et al.l [l993). A com- 
plete list of pu blication s from this project can be found at 
http://www.usm.uni-muenchen.de/people/ip/papers/papers. 
Within this project we have been systematically working 
on the data for the low ionization stages of iron peak 



Table 1. Configuration expansion for Can, and scal- 
ing parameters X n i for each orbital of the configurations 
ls 2 2s 2 2p 6 3s 2 3p 6 nZ in the Thomas-Fermi-Dirac potential 

nl Configurations 



3d. 4s, 4p, 4d, 4f, 
6g, 7s, 7d, 7f, 7g, 



5s, 5p, 5d, 5f, 5g, 6s, 6p, 6d, 6f, 
8s, 8d, 8f, 8g 



elements, 



e^ 



radiative and collisional rate s for Fe i-iv 

20041) . 



Bautista & Pradhanl Il998f). Niii dBautistal 
Bautistal l200ll ), Niiv dMelendez fe Bautista 



Ni in 
I2005D . The 

objective of the present work is to provide accurate and 
complete atomic data for a detailed spectral modeling of 
Call. Such a model should be large enough to account 
for various processes such as collisional excitation includ- 
ing cascades from high levels, fluorescence by line and 
continuum radiation, and line optical depth effects. 



, Is:, 1.43880, 2s: 1. 
ypT 1.09850, 3d: 1 
4d: 1.07690, 4f: 1. 
5d: 1.07660, 5f: 1. 
6p: 1.07820, 6d: 1 
7s: 1.08510, 7d: 1. 
8s: 1.08590, 8d: 1. 



11310, 2p: 
.07950, 4s: 
05000, 5s: 
04950, 5g: 
.07690, 6f: 
.07760, 7f: 
.07860, 8f: 



1.05670, 3s: 
1.08770, 4p: 
1.08510, 5p: 
1.01940, 6s: 
1.04960, 6g: 
1.05010, 7g: 
1.05120, 8g: 



1.10580, 
1.07730, 
1.07690, 
1.08480, 
1.01920, 
1.01930, 
1.01950 



Table 2. Term energies for Call (in Ryd). The table shows 
results computed without PI (w/o PI), with PI and exper- 
imental energies from NIST V.3.1.0 



2. Atomic data 

2.1. Atomic structure calculations 

We use t he at o mic st ructure code AUTOSTRUCTURE 
(|Badnelll Il986l Il997t ) to reproduce the structure of 
the Can ion. This code is based on the pro- 
gram SUPERSTRU CTURE originally developed by 
lEissner et al.l (|1974| ). but incorporates various improve- 
ments and new capabilities like two-body non-fine-structure 
operators of the Breit-Pauli Hamiltonian and polarization 
model potentials. In this approach, the wave functions are 
written as configuration interaction expansions of the type: 

if>i = ^2<f>jCji, (1) 
3 

where the coefficients Cji are determined by diagonaliza- 
tion of I H I ipi). Here H is the Hamiltonian and the 
basic functions <pj are constructed from one-electron or- 
bitals g enerated using the Thomas- Fermi-Dirac model po- 
tential (|Eissner fe Nussbaumer|[l969l ) . including A„; scaling 
parameters which are optimized by minimizing a weighted 
sum of energies. The basic list of configurations and scaling 
parameters used in this work are listed in Table 1. 

Relativistic effects are included in the calculation by 
means of the Breit-Pauli operators in the form: 

H = H nx + i?bp, (2) 

where H nr is the usual non-relativistic Hamiltonian and 
-ffbp is the Breit-Pauli p erturbation, which includes one- 
and two-body operators (|Jonesl I1970L 119711 : lEissner et al.l 
119741) . 

2.2. Model potential 

In order to obtain accurate orbitals in our multiconfigura- 
tion frozen-core approximation we include the polarization 
interaction between the valence electron and the core in a 
model potent ial. We used a mode l poten tial V po i of the form 
described by iNorcross fe Seatonl (Il976f) : 

V pol {r lP ) = -^[l-eM-r/pf], (3) 



TERM w/o PI PI NIST 



where ay is the static dipole core polarizability of the ion 
Cain and p is adjusted empi rically to yield good agreement 
with experimental energies. IWallerl {l926) using the "non- 
penetrating" orbitals theory obtained ad = 3.31 for the d, 
f and g states. We adopt p = 2.25 the cut-off parameter. 
This yields accurate binding energies for the n=3,4,5,6,7 
and 8 configurations of Can. 

The expansion considered here for the Can system in- 
cludes 23 LS terms. Table 2 presents the complete list of 
states included as well as a comparison between the calcu- 
lated and observed target term energies, averaged over fine 
structure. Here, we show the energies without polarization 
interaction (w/o PI) and those with polarization interac- 



1 


4s 


2 S 


0.000000 


0.000000 


0.000000 


2 


3d 


2 D 


0.147449 


0.124596 


0.124721 


3 


4p 


2po 


0.219659 


0.232385 


0.230916 


4 


5s 


2 S 


0.454729 


0.478864 


0.475380 


5 


4d 


2 D 


0.500830 


0.524648 


0.518062 


6 


5p 


2po 


0.529136 


0.555242 


0.552092 


7 


4f 


2po 


0.593486 


0.623006 


0.620180 


8 


6s 


2 S 


0.618814 


0.647386 


0.644061 


9 


5d 


2 D 


0.639097 


0.667632 


0.662741 


10 


6p 


2po 


0.652978 


0.682121 


0.678980 


11 


5f 


2po 


0.683576 


0.713962 


0.711101 


12 


5g 


2 G 


0.683935 


0.715180 


0.712289 


13 


7s 


2 S 


0.697115 


0.727147 


0.723986 


14 


6d 


2 D 


0.707827 


0.737833 


0.733791 


15 


6f 


2po 


0.732574 


0.763404 


0.760526 


16 


6g 


2 G 


0.732824 


0.764166 


0.761272 


17 


8s 


2 S 


0.740613 


0.771273 


0.768206 


18 


7d 


2 D 


0.746957 


0.777594 


0.773987 


19 


7f 


2po 


0.762129 


0.793202 


0.790315 


20 


7g 


2 G 


0.762303 


0.793703 


0.790807 


21 


8d 


2 D 


0.771345 


0.802302 


0.798935 


22 


8f 


2po 


0.781312 


0.812527 


0.809636 


23 


8g 


2 G 


0.781436 


0.812872 


0.809975 
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Fig. 1. log gfv plotted against log <?/l for transitions be- 
tween energy levels. Panels (a) and (b) show the results 
computed without polarization interaction and with polar- 
ization interaction respectively 



tion (PI). It can be seen that the contribution of PI can 
reach up to 15%, especially for the lower energy terms. 

In the calculation of radiative rates, fine tuning of eigen- 
states is performed with term energy corrections (TEC), 
where the improved relativistic wave function, i/jf 1 , is ob- 
tained in terms of the non-rclativistic functions 



LS 



| H bp | #*) 



E 3 S 



(4) 



with the LS energy differences E^ s — E^ s adjusted to fit 
weighted averaged e nergies of the experimental multiplets 
(jZeippen et al.lll977ft . 

In our best target representation, which accounts for the 
interaction between the valence electron and the core, the 
theoretical energies for all the 23 terms are typically within 
2% of the experimental values before any further empirical 
correction. After TEC, the agreement with experimental 
energies is better than 1%. 

For dipole-allowed transitions, spontaneous decay rates 
are given by 



1 



Af = 2.6774 x 10 9 (^ - E.f-Sf/ (s' 1 ), 



fj, 



(5) 



while for forbidden transitions we consider electric 
quadrupole (E2) and magnetic dipole (Ml) transition rates 
given by 



1 



Af 2 = 2.6733 x 10 3 (K t - Ejf-S? 2 (s' 1 ) 



and 



A™ 1 = 3.5644 x 10 4 (£, 



yi 



(6) 



(7) 



Here, gi is the statistical weight of the upper initial level i, 
Sij is the line strength and E is the energy in Rydbergs. 

Eqns. (|5l6l and[7| show that the transition rates are sen- 
sitive to the accuracy of the energy levels, particularly for 



forbidden transitions among nearby levels. Thus, we per- 
form further adjustments to the transitions rates by cor- 
recting our best calculated energies to experimental values. 

In Fig. [1] we plot the gf- values for dipole allowed tran- 
sitions among fine structure levels computed in the length 
gauge vs. those in the velocity gauge. We present the gf- 
values without PI (a) and with PI (b). The overall agree- 
ment between the two gauges is around 5% for log(g/)- 
values greater than —3 when accounting for PI and greater 
than 15% without PI. This is a good indicator of the quality 
of the dipole allowed radiative data. 

In Table [3] we present an extensive comparison between 
the present results and previous lifetimes for the metastable 
levels 3d 2 D 5 / 2 and 3d 2 D 3 / 2 . These levels are of particu- 
lar astrophysical interest because they yield the prominent 
spectral lines AA 7293, 7326 A. Our results including PI and 
TEC are in excellent agreement with experimental values, 
while the results that neglect PI are ~10% too low. 

Table 3. Lifetimes of the metastable 3d 2 D levels of Can 



Level 



Present 
w/o PI PI 



Other Experiment(r(s)) 



3d 2 D 3/2 0.926 



3d 2 D 5/2 0.901 



1.107 



1.08 



1.081 1 

1.16 2 

1.27 3 

0.98 4 

1.196 5 

1.058 1 

1.14 2 

1.24 3 

0.95 4 

1.165 5 

1.045 6 

1.14 11 



1.176±0.011 5 



1.168±0.009 5 

1.152±0.020 7 

1.100±0.018 8 

1.054±0.061 9 

1.149±0.014 10 

1.064±0.017 n 



^Zeinpenl (l!990h. Waeck et all (Il992h. 3flnet fc Johnson! 
(119911) . lAli feKiml (11988IL iKreuter et all (l2005lL 
^Liawi (11995T). iKnoon et all (l2004h. iBlock et all 
111999ft. HArbes et aU (|1994D . lc (Staanum et all (l2004h . 
1] lGudions et alJ (| 1996ft 



In Table[4]we compare the calculated lifetimes for short- 
lived levels of Ca II from the present calculations with other 
theoretical and experimental values. For the lower levels 
[Ap 2 P°/ 2 an d 4p 2 -P3/ 2 ) the effect of polarization interaction 
is ^20%. Overall, the differences between the results of our 
best model and experimental values are less than 5%, ex- 
cept for the level 5d 2 D 5 / 2 - For this level, our result agrees 
with previous independent calc ulations but is about 4 0% 
below the experimental values of lAndersen et al"] (|1970ft . A 
new measurement of this lifetime would be very important. 

In Table [5] we present a comparison between calcu- 
lated and experimental oscillator strengths in absorption. 
For the calculated oscillator strengths we choose the most 
complete and representative set of values as well as the 
most cited. We c ompare with previous calculations by 
IVaeck et all (|1992ft where they use a multiconfiguration 
Hartree-Fock method with core polarization included vari- 
ationally (SECP) using a model potential of the form de- 



4 



M. Melendez et al.: Atomic data from the Iron Project LXIV 



Table 4. Can lifetimes (in ns). The second column gives the results with neither PI (w/o PI) nor TECs (w/o TECs), 
the third column gives the results without PI (w/o PI) but with TEC, the fourth columns gives results with PI but no 
TEC, and the fifth column shows the results with both PI and TEC. Theoretical (Other) and experimental (Experiment) 
values from other authors are given in the sixth and seventh columns, respectively. 



Level Present 
w/o PI 
w/o TEC TEC 



PI 

w/o TEC TEC 



Other Experiment 



4p 


2po 
r l/2 


6.978 


5.734 


6.697 


6.837 


6.44° 


7.07±0.07 b 












6.87 b 

6.39 c 

6.94 d 

7.045 e 

7.047' 


7.5±0.5 S 
6.62±0.35 h 
6.9^0.^ 
7.098±0.020 


4p 


2po 
r 3/2 


6.797 


5.577 


6.508 


6.649 


6.28° 


6.87±0.06 h 












6.24 b 

6.21 c 

6.75 d 

6.852 e 

6.833' 


7.4±0.6 9 

6.68±0.35 /l 

6.72±0.20 J 

6.61±0.30 fc 

6.87±0.18 i 

6.924±0.019 


4d 


2 D 3/2 


2.963 


2.779 


2.781 


2.934 


2.868 e 




4d 


2 D 5/2 


2.981 


2.799 


2.800 


2.952 


2.886 e 




4f 


2po 

5/2 


3.625 


2.656 


3.487 


3.451 


3.895 e 




4f 


2po 
r 7/2 


3.625 


2.654 


3.486 


3.448 


3.897 e 




5s 


2 ^l/2 


4.310 


3.833 


3.886 


3.982 


4.13° 


4.3±0.49 












4.153 e 
3M f 




5p 


2po 
r l/2 


40.195 


34.903 


36.541 


35.174 


33.78° 














36.200 6 




5p 


2po 
r 3/2 


38.987 


33.755 


35.783 


34.401 


33.92 a 














35.349 e 




5d 


2D 3/2 


5.960 


5.625 


5.895 


6.102 


6.148 e 




5d 


2 2 D 5/2 


6.002 


5.669 


5.944 


6.154 


6.199 e 


4.3±0.2f 


Gs 


Sl/2 


7.007 


6.390 


6.391 


6.457 


6.90 a 














6.766 e 
6.51-f 




6p 


2po 
r l/2 


114.06 


98.569 


93.141 


90.173 


58.95 a 














100.254 e 




6p 


2po 
r 3/2 


111.23 


95.788 


92.262 


89.129 


43.06 a 
99.675 e 





a Hafncr & Schwarz 


(119781). f Gosselin et al.l (Il988a). 'Wiese et al. 


(1969). iGuet & Johnson 


(19911). iTheodosiou (|1989h 


* Black et al. (1972 

k-ic\cc& tn.j o 


). ^Andersen et all (1970 

_ o„i 1 \-t cvich At : J / 


). lAnsbacher et al. 

1 nnch ttA t: P_ /"ii 


(1985). iGosselin et al.l ( 


1988H). 4Smith & Gallaeherl 



scribe d first bv lBavlisl l|l977t ) and then reviewed b y lHibbertl 
(|l989h . We also compare with iTheodosioul (|l989l ) who used 
a Hartree-S later (HS) core potent ial with semiempirical 
corrections, iGuet fe Johnsonf (ll99ll) who used relativistic 
many-body theory and lBrage et al.l (| 19931 ) using an orthog- 
onal Breit-Pauli MC HF method. We also compare with ex- 
perimental values of lGallagherl (|l967t ) who used the Hanle- 
effect technique with optical excitation from the ion ground 
state. The inclusion of PI can affect the oscillator strengths 
by up to 12%. One finds very good agreement between the 
present calculation, previous calculations and experimental 
values. 



In Table [5] we present our electric quadrupole and mag- 
netic dipole transitions probabilities computed including 
P I and TEC. We compare these with previous calculations 
of IZeippenl (|1990[ ). who used a two-step minimization pro- 
cedure and semi-empirical term energy correc tions using 
the co mputer program SUPERSTRUCTURE, lAli fc Kiml 
(|1988[ ) who used the multic onfiguration Dirac-Fock formal- 
ism and IVaeck et al.l (|1992l). Our r esults are in good agree- 
ment (within ~2%) wit h IZeippenl (1 19901) a nd within ~13% 
and ~ 8% with respect to lAli fc Kiml (|l988t ) and lVaeck et alj 
((1992T) . 
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Table 5. Oscillat o r strengths for Can. T he o scillator strengths from previous calculations of Vaec k et al.l (fl 992). 
iTheodosioul (|1989h . iGuet fe; Johnson! (|1991h and iBrage et alj (|1993h . and experimental data from iGallagherl (|1967l ) are 
also given. 

Transition Present Vaeck Theodosiou Guet Brage Experiments 
w/o PI PI 

4s Si/2 - 4p P° /2 0.364 0.323 0.318 0.316 0.320 0.321 



4s 2 S 1/2 - 4p 2 P° /2 0.734 0.652 0.641 0.637 0.645 0.649 0.66±0.02 

3d 2 D 3/2 - 4p 2 P° /2 0.0478 0.0537 0.0547 0.0473 0.0494 0.0524 

3d 2 D 3/2 - 4p 2 P° 0.0098 0.0110 0.0112 0.0096 0.0101 0.0107 0.0088±0.001 



3/2 
1/2 

J 3/2 - ± V -r>° /2 

) 5/2 - 4p 2 P° /2 0.0584 0.0656 0.0666 0.0574 0.0601 0.0637 0.053±0.006 



Table 6. Electric quadrupole (E2 ) and magnetic dipole (Ml) t ransit ions probabilities of [Ca ll], We compare our results 
with previous calculated values of lZeippenl (|1990D . lAli fc Kiml (|1988T ) and ! Vaeck et all (|1992D . 



M s ~ L ) 




Type 


Present 


Zeippen 


Ali and Kim 


Vaeck 


3d 2 D 3/2 
3d 2 D 5/2 
3d 2 D 5/2 


- 4s 2 S 1/2 

- 4s 2 S 1/2 

- 3d 2 D 3/2 


E2 
E2 
Ml 


0.905 
0.928 

2.41xl0" 6 


0.925 
0.945 

2.45x10^ 


1.02 
1.05 

2.45xl0~ 6 


0.84 
0.86 



2.3. Scattering calculations 

In the close coupling (CC) approximation the total wave 
function of the electron-ion system is represented as 

iP(E;LStt)=AJ2xA + J2 c ^ ( 8 ) 

i 3 

where \i is the target ion wave function in a specific state 
LSi, 0i is the wave function of the free electron, are short 
range correlation functions for the bound (e+ion) system, 
and A is the antisymetrization operator. 

The variational procedure gives rise to a set of coupled 
integro-differentia l equations that are solved with the R- 
matrix technique (jBurke et al.ll97ll : lBerrington et al.ll978t 
Il995h within a box of radius r < a. In the asymptotic re- 
gion, r > a, exchange between the outer electron and the 
target ion can be neglected and if all long-range poten- 
tials beyond Coulombic are also neglected, the reactance K- 
matrix and the scattering S'-matrix are obtained by match- 
ing at the boundary the inner-radial functions to linear 
combinations of the outer-region Coulomb solutions. Later, 
the contributions of long-range poten tials to the collisio n 
strengths are included perturbatively (jGriffin et al.| [T998). 

We use the LS'-coupling R-matrix method that in- 
cludes mass-velocity and Darwin operators and the po- 
larization model potential. Note that the scattering cal- 
culations include one-body relativistic operators while the 
atomic structure calculations using AUTOSTRUCTURE 
include two-body operators of the Brcit-Pauli Hamiltonian 
namely (two-body) fine-structure operators and (two-body) 
non-fine-structure operators. 

The ^-matrix elements determine the collision strength 
for a transition from an initial target state i to a final target 
state /, 

where w = (2L + 1)(2S + 1) or (2J + 1) depending on the 
coupling scheme, and the summation runs over the partial 



waves and channels coupling the initial and final states of 
interest. 

In order to derive fine-structure results out of the LS- 
coupling calculation we employ the intermed iate-coupling 
frame transformation (ICFT) method of iGriffin et alj 
(jl998l ). The ICFT method uses the multi-channel quantum 
defect theory (MQDT) to generate the LS-coupled unphys- 
ical JC-matrices. In this approach one treats all scattering 
channels as open and calculates the term-coupling coeffi- 
cients (TCCs) to transform the unphysical -ftf-matrices to 
full intermediate coupling. Finally, we can generate the 
physical fC-matrices on a fine energy mesh. Because all 
channels are treated as open this method eliminates the 
problems associated with the transformation of the physi- 
cal S'-matrices with closed channels and consequently yields 
accurate results for both background and resonances of the 
collision strengths at all energies. 

The computat ions were carried out w ith the RMATRX 
package of codes ()Berrington et al1l!995h and also include 
the dipole polarization potential for the interaction of the 
valence electron and the core. The set of (N+l)-electron 
wave functions to the right of the CC expansion in Eq. (jHJ) 
includes all the configurations that result from adding an 
additional electron to the target configurations. 

In Fig. [5] we compare the collision strengths for differ- 
ent target expansions. For the various close-coupling expan- 
sions we include all the terms that arise from the n=3, 4, 5, 
6, 7 and 8 configurations. It is interesting to note that by 
increasing the number of terms in the close-coupling expan- 
sion the resonances tend to become better organized and 
blended, leading to regular broad series of structures. This 
is because of the increasing number of channels for decay 
of autoionizing levels. We adopt the expansion up to n=8 
for all further calculations. 

Partial wave contributions to the summation in Eqn. 3 
were included from 162 SLir total symmetries with angular 
momentum L = — 40, total multiplicities (25 + 1) = 1 
and 3, and parities even and odd. As for the number of 
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Fig. 2. Comparison between collision strengths for 
the Can ion with different CC expansions for the 
3d 2 D 3/2 - 4s 2 S 1/2 (upper panel) and 4p 2 P° /2 - 4s 2 S 1/2 
(lower panel) transitions, n represents the principal quan- 
tum number for the valence electron. 



continuum orbitals per angular momentum we found nec- 
essary to include up to 40. The collision strengths were 
"topped up" with estimates of the contributions of higher 
partial waves for the optically allow ed trans i tions based on 
the Coulomb-Bethe approximation Bur gess! (|1974f ) and for 
non-allowed transitions we approximate the top-up with a 
geometric series. It was verified that by explicitly including 
partial waves up to L = 40 such "top ups" amount to less 
than 10% of the total collision strengths for all transitions 
with the exception of the transitions ng-ni (for n >6) for 
which the top up is greater than 20%. This is because of 
the very small energy difference between these terms. 

In order to check for convergence of the partial wave ex- 
pansion one needs to study the high energy behavior for the 
collision strengths (f2). In Fig. [3] we plot reduced collision 
strengths (f2 r ) as function o f reduced energy (E r ) f ollow- 
ing the procedure describe bv lBurgess & Tullv (1992). This 
approach allows us to visualize the complete range of ener- 
gies mapped onto the interval [0,1]- The last point in these 
graphs represents the infinite energy limit of the collision 
strengths (E oo). For this plot we adopt a scaling pa- 
rameter C = 1.5. These plots show good progression of the 
collision strengths toward the high energy limit, which gives 
confidence on the consistency and quality of the data. 

Fig. 2] shows the collision strength computed with and 
without polarization interaction for a sample of transitions 
between levels of the lowest four multiplets. One can see 
that the core polarization has an important effect on the 
collision strengths. We verified through various calculations 
that such effects come almost entirely from the change in 
the radial wave functions of the target, while the effects of 
the dipole polarization operator in the scattering matrix 
is nearly neg l igible . The same conclusion was derived by 
iMitrov et aU(|l988D . 

The collision strengths were calculated at 22000 energy 
points from to 4 Ry, with a resolution of 10~ 5 Ry in the 
region with resonances and 5 x 10 -3 Ry at higher energies. 
This number of points was found sufficient to resolve most 



Fig. 3. Reduced collision strengths for the 4s 2 S — 4p 2 P° 
(upper panel), 3d 2 D — 4p 2 P° (middle panel) and 
4s 2 S — 3d 2 D (lower panel) transitions. 



resonance structures for accurate calculations of effective 
collision strengths for temperatures above ~5000 K. 

A dimensionless thermally-averaged effective collision 
strength results from integrating the collision strength over 
a Maxwellian distribution of electron velocities 



n i/ exp(-e / /fcT)d(e//fcT), 



(10) 



where e/ is the kinetic energy of the outgoing electron, T 
the electron temperature in Kelvin and k = 6.339 x 10~ 6 
Ry/K is the Boltzmann constant. 

In Table [7] we compare the present effective collision 
strengths for the 4s — > 4p transition of Call. W e show 
the values obtained bv lOsterbrock fc Wallace! (119771) wher e 
they used experimental results of iTavlor fc DunrT (ll973l), 
as reduced to collision str engths bvlSeatonl (1 19751) . We also 
present the values from IBurgess et al.l ( 1995|) as we use 
their five point cubic spline fitting param eters to tabulate 
the T (T) using the procedure described in lBurgess fc Tullvl 
(1992). We compare them with our present calculation 
without PI in the target representation and in the scat- 
tering calculations (w/o PI); without PI in the R-matrix 
calculations o nly (PI-S) and with PI. Our b est results agree 
within 11% of lOsterbrock fc Wallace] (Il977l) valu es. On the 
other hand, the results of IBurgess et al.l (|l995l ) obtained 
with the non-exchange distorted wave approximation ap- 
pear overestimated by ~50%. In this calculation, the Call 
target is represented in the frozen core approximation ne- 
glecting polarization. As we pointed out before the main 
contribution of the dipole polarization potential is in the 
representation of the target ion. Further, our effective col- 
lision strengths when neglecting polarizati on are overes- 
timate d by 30%, in closer agreement with IBurgess et al.l 
(fl995h . 



3. Conclusions 

We have computed radiative data, collision strengths and 
effective collision strengths for transitions among 41 levels 
from the n=3, 4, 5, 6, 7 and 8 configurations of Call. The 
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Fig. 4. Comparison between collision strengths for the Call ion. The solid and dashed lines depict the results with and 
without model potential interaction respectively. The transitions are: (a) 3d 2 D 3 / 2 — 4s 2 Si/ 2 ; (b) 3d 2 D 5 / 2 — 4s 2 Si/ 2 ; 
(c) 4p 2 P° - 4s 2 S 1/2 ; (d) 4p 2 P° - 4s 2 S 1/2 ; (e) 3d 2 D 3/2 - 4p 2 P° and (f) 3d 2 D 5/2 - 4p 2 P° . 



Table 7. Effective collision strengths for the 4s 2 S - 4p 2 P° 
transition Can. The first three lines correspond to results 
of the present calculations for the cases of no PI in cither 
the target orbitals of scattering calculation (w/o PI), PI in 
the scattering calculation only (PI-S), and PI in the tar- 
get and scattering calculations (PI) . These results are com- 
pared with values deduced fro m experimental cross sections 
ijOsterbrock & Wallace] I1977D COW") and th e most recent 
theoretical calculations ([Burgess et aLlU^ S) (BCT). 





T(k) 






5000 


10000 


15000 


20000 




19.03 


20.82 


22.45 


24.03 


w/o PI 


16.96 


19.34 


21.35 


23.13 


PI-S 


17.07 


19.44 


21.45 


23.23 


PI 


15.6 


17.5 


19.2 


20.8 


OW 


24.87 


27.50 


29.88 


32.04 


BCT 



radiative data were calculated using the Thomas-Fermi- 
Dirac central potential with a model core potential that 
account for dipole polarization interaction of the valence 



electron with the core. We also present an extensive com- 
parison between our results and the most recent experi- 
ments and calculations for the lifetimes of Can. 

Effective collision strengths are available for various 
temperatures that expand from 3000 K to 38000 K. 
The whole set of data reported here including en- 
ergy levels, infinite energy limit Born collision strengths, 
transition probabilities and effective collision strengths 
can be obtained in electronic form at the CDS via 
anonymo us ftp to cdsarc. u-strasbg.fr (130.79.128.5), via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat7J/A-l-A/ or by re- 
quest T!o - fJ^e - atitliOTs! 
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